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Abstract 



There are recent interests with CsI(Tl) scintihating crystals for Dark Matter 
experiments. The scattering signatures by neutrons on a CsI(Tl) detector were 
studied using a neutron beam generated by a 13 MV Tandem accelerator. The 
energy spectra of nuclear recoils from 7 keV to 132 keV were measured, and their 
quenching factors for scintillating light yield were derived. The data represents the 
first confirmation of the Optical Model predictions on neutron elastic scatterings 
with a direct measurement on the nuclear recoils of heavy nuclei. The pulse shape 
discrimination techniques to differentiate nuclear recoils from 7-background were 
studied. Internal consistencies were obtained among the different methods of light 
yield measurements. The projected capabilities for Cold Dark Matter searches with 
CsI(Tl) crystals are presented. 

PACS Codes: 25.40.Dn, 95.35. +d, 29.40.Mc. 
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1 Introduction 



The detection of Dark Matter and the studies of their properties [|l| are of fundamental 
importance in particle physics and cosmology. The Weakly Interacting Massive Particles 
(WIMPs) are good candidates for "Cold" Dark Matter, and their experimental searches 
have gathered a lot of interests in recent years. The most promising avenue is to detect 
the nuclear recoil signatures due to elastic scatterings of WIMPs on the target isotopes. 
The typical energy depositions are only of the order of 10 keV, imposing big experimental 
challenges in terms of the detection of weak signals as well as background control at low 
energy close to detection threshold. A wide spectrum of experimental techniques is being 
pursued [^. 

The DAMA experiment observed an annual modulation of nuclear recoil events 0] 
with Nal(Tl) scintillating crystal detectors, which can be interpreted as positive evidence 
of WIMPs due to the difference of the relative velocities of the Earth from the Halo sea 
within the year. This, however, contradicts the limit from the CDMS experiment based 
on cryogenic technique. 

There is still much room for new detector concept to thoroughly test the DAMA 
parameter space, and to push the sensitivities further. It would be of great interest if 
the sensitivities of WIMP searches can probe the level predicted by the various Super- 
Symmetry models. 

There are potential merits of using CsI(Tl) scintillating crystals for WIMP search 
and other low-energy low-background experiments ||^, |^. An experiment towards 200 kg 
of CsI(Tl) crystal scintillators to study low energy neutrino interactions at the Kuo-Sheng 
power reactor is being pursued [[7|, while the adaptation of the crystal for Dark Matter 
searches are the focus of several on-going projects p], H, p!0| . 

The high-A content of the Csl enhances the sensitivities for the spin- independent in- 
teractions (which depends on A^) between the WIMPs and the target, relative to most 
other candidate target isotopes. The high-Z composition allows a compact design and 
provides large suppression of background due to ambient radioactivity if a three dimen- 
sional fiducial volume definition can be realized. Both ^^^Cs and ^^^I are 100% in their 
respective isotopic abundance. Being close in their mass numbers, the response to nuclear 
recoil from the interactions with WIMPs would be similar, allowing simpler interpretation 
of experimental signatures. 

As a detector, the crystal has large light yield, low energy threshold and with pulse 
shape discrimination characteristics for 7/0; separation. Its characteristics of being only 
slightly hygroscopic implies that it can be machined easily and does not require hermetic 
seal (that is, passive materials) in a large detector system. In addition, large (40 tons) 



electromagnetic calorimeter systems [H have been constructed and made operational in 
high energy physics experiments, making this technology affordable and realistic to scale 
up. Considering all the associated costs, the price of CsI(Tl) is in fact less than that for 
Nal(Tl). In order to produce positive and definite evidence of the WIMPs, an accurate 
measurement of the annual modulation (where the maximal effects are only 7%) would 
be necessary such that the availability of large target mass is a very desirable feature. 

One of the key issues to realize a Dark Matter search experiment with CsI(Tl) crystal 
scintillator is the studies of the experimental signatures of nuclear recoils due to WIMP- 
nuclei elastic scatterings. Nuclear recoils produce large charge density (dE/dx) such that 
the scintillating light yield is "quenched" and the timing profile of pulse shape is different 
relative to the same energy deposition by minimum ionizing particles [l^ . 

These signatures are the same as the nuclear recoil events produced by elastic scat- 
tering of neutrons on nuclei. A measurement was performed with the neutron facility at 
the 13 MV Tandem accelerator at the China Institute of Atomic Energy at Beijing. The 
results of the quenching factors measurement and of pulse shape discrimination studies 
are reported in this article. They extend the work from the other recent measurements on 
CsI(Tl) H], as well as on other scintillating crystals such as Nal(Tl) and CaF2(Eu) ^ 



2 Experimental Setup and Procedures 

The experiment was performed at HI- 13 tandem accelerator at the China Institute of 
Atomic Energy (CIAE) in Beijing. A schematic diagram is shown in Figure |^. A pulsed 
deuteron beam at 5.6 MeV interacted with a deuterium gas target in a cell at a pressure 
of 6 bar and the dimensions of 1 cm in diameter and 3 cm in length. Neutrons at 8 MeV 
kinetic energy were produced at zero degree. The CsI(Tl) sample was located 2.02 meters 
away from the deuterium target. Neutrons at zero degree were selected by a 32 mm by 
35 mm collimator of length 1.2 m. The collimator were surrounded by shielding materials 
like iron, paraffin, lead and polyethylene to reduce background. During this data taking, 
the repetition rate of the pulsed beam was set at 4 MHz and the average beam current 
was at 1 yuA. 

The scattering target, which also functioned as a detector, was a CsI(Tl) crystal 
scintillatoiQ with diameter 3 cm and length 3 cm wrapped with teflon sheets, aluminium 
foil and black plastic tape. To minimize supporting structures, the detector was hung at 
the correct position by a piece of string. The readout was achieved by a 29 mm diameter 
photo-multiplier tube (PMT)| The PMT signals passed through an amplifier and shaper 

'''Producer: Unique Crystal, Beijing 
■'■CRllO, Hamamatsu Photonics, China 



and were digitized by a 20 MHz, 8-bit, Flash Analog-to-Digital Convertor (FADC) [0 
developed for the Kuo-Sheng reactor neutrino experiment 0. The trigger was provided 
from a time of flight (TOP) system to tag the elastically scattered neutrons. The 7- 
response of the CsI(Tl) detector was calibrated with an LED pulser operated at single 
photo-electron intensity, as well as standard ^°^Cd and ^^'^Ba 7-sources. The calibration 
is 4 photo-electrons per keV of electron-equivalence energy. 

Neutron tag was provided by liquid scintillatoi]^ detectors 105 mm in diameter and 
50 mm in length equipped with PMTQ readout. The liquid scintillator provides pulse 
shape discrimination (PSD) capabilities for n-7 differentiation. The PSD was achieved by 
commercial electronics£[ The neutron tag signal was used to define the "START" timing 
for the TOP system. The pulsed deuteron "pick-off" signal provided by the accelerator 
was delayed and used as "STOP" for the TOP. The neutron detectors were placed at a 
distance of 0.5 to 2 m from the CsI(Tl) target at different angles, to provide optimized 
TOP resolution to differentiate events due to neutron elastic scattering off hydrogen atoms 
or heavier nuclei in the wrapping materials. 

The PADC recorded data continuously in a circular buffer of 4k size. A START-STOP 
sequence from the TOP system initiated a trigger which stopped the PADC digitization 
after 25 /is. The pulse shape from the CsI(Tl) crystal was recorded for a pre-trigger and 
post-trigger periods of 5 /is and 25 /is, respectively. With the timing settings at data 
taking, the CsI(Tl) pulse started at —1.25 /is if it was due to nuclear recoils of heavy 
nuclei. This timing information provides a powerful means of the rejection of background 
due to proton recoils and pile-up events. 

Data taking at each scattering angle was complemented by a measurement of the 
TOP background spectrum with empty target which was made of the CsI(Tl) wrapping 
materials only. Displayed in Pigure ^ is the TOP spectra at 50° with empty and complete 
target. The lower and upper peaks are due to neutron elastic scattering off protons and 
the heavier nuclei such as carbon, aluminium, cesium, and iodine, respectively. The fiat 
background is due to random coincidence. It can be seen that with the CsI(Tl) target in 
place, the upper TOP peak for neutron scattering off heavier nuclei becomes enhanced. 
The difference between the two TOP spectra is the nuclear recoils spectra from Cs and I. 
The proton-recoil events tagged by the low TOP peak provides a means to evaluate the 
signal yield and to confirm the distance between neutron detector and Csl target. 

The reconstructed energy is defined as the pedestal-subtracted integrated area of the 
PADC pulse over 2.5 /is starting from the to of PADC (1.25 /is before TOP trigger in 
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this data taking). This restricted time window helps to minimize the pile up effect from 
accidental low energy 7's which are abundant along the beam line. Events are categorized 
into two groups by the central value between the two TOF peaks and their corresponding 
nuclear recoil energy spectra are compared. The energy spectra at 50° are shown in 
Figure ^ and for events with high TOF values and low TOF values, respectively. 

The energy spectrum tagged by the low TOF trigger is displayed in Figure Since 
the trigger is due to proton recoils from the wrapping materials, there should be no corre- 
lated signals from the CsI(Tl) target. The energy spectrum which shows a distinguished 
peak with long tail on the high-energy side is due to events from accidental 7's arriving 
within the integration time-window. 

The detector response for Csl nuclear recoil events was obtained by subtracting Fig- 
ure from Figure |^ with proper normalizations, where the factors are derived by nor- 
malizing the low TOF peaks from proton recoils between the full and empty target spectra. 



3 Measurement Results 

The resulting Csl recoil spectra obtained by the "statistical background subtraction" 
scheme discussed in Section ^ are shown in Figure ^ for eight different scattering angles. 
It can be seen that this procedure is valid since the background levels away from the peaks 
after subtraction are consistent with zero. 

The nuclear recoil energy (T) is related to the neutron scattering angle [9) by a simple 
kinematical formula: 

2A 

T = (1 - cos^) E„ , (1) 

where A is the target atomic mass and En is the incident neutron energy. The neutron 
scattering angle can be considered to be the same in this case in both the laboratory 
and center-of-mass systems. Comparisons of the nuclear recoil spectra with those due to 
calibration sources allow the evaluation of the quenching factors, defined as the light yield 
from nuclear recoils versus that from 7 sources. 

To enhance the signal-to-noise ratio in the large background environment, the nuclear 
recoil pulses are integrated for a restricted 2.5 fis (or 50 FADC time bins), as compared 
to the full time window 12.5 /is (or 250 FADC time bins) for 100% light collection in 
the 7-sources calibration measurements. A correction factor of 1.15 is applied to the 
nuclear recoils light yield to account for this partial summation. This factor was derived 
by averaging and comparing a large number of events due to nuclear recoils and ^°^Cd 
7-source, as depicted in Figure The estimated uncertainty is 10%, based on studies 
with different integration time window from 1.5 /is to 12.5 /is. 



The deviation from the central value of nuclear recoil energy due to finite detector size 
is at most 2% for different scattering angles. This is checked by calculating the average 
recoil energy with proper weights of event angular distribution and detector acceptance. 
A shift of detector at 1 mm scale has negligible effect on the central value of nuclear recoil 
energy. The beam current was stable at the 1% level during data taking, as indicated 
by the stable data taking rate (typical value is about 4 Hz) as well as readings from 
beam monitor counters. The main systematic uncertainties originate from the analysis 
procedures, and amount to 16% at the lowest energy point at 20°. The statistical and sys- 
tematic uncertainties are then combined in quadrature to provide the total uncertainties 
for each data point. 

The results at different scattering angles are summarized in Table |[ The measured 
quenching factors are depicted in Figure while the variation of electron-equivalence 
light yield (the "visible" energy) versus recoil energies are shown in Figure Ipo. The 
other measurements P, P] are overlaid for comparison. There is a clear trend of less 
quenching towards low energy points. This measurement achieved a lower threshold 
and with improved uncertainty compared to Ref. 0. The quenching factors are of the 
similar range to those for Nal(Tl) [|15|, which are typically 25% for Na and 8% for I. 
It can be seen that the light yield is higher for interactions with the high-A isotopes 
in CsI(Tl), providing big advantages to probing the spin-independent interactions. In 
addition. Figure |]b depicts a linear relationship between the electron-equivalence light 
yield and the nuclear recoil energy, as indicated by a linear fit our data set. Data from 
Ref. shows more quenching and deviates from the linear regime at high energy. This 
is different from the highly non-linear response due to threshold effects observed in liquid 
scintillator [ffH] at the same range of recoil energy. 

The total number of observed signal events are given by the area under the Gaussian 
fits in Figure ^ To derive the differential cross sections, the data taking time, the ac- 
ceptance (proportional to r~^) and the efficiency correction factors have to be evaluated. 
The efficiency factors are due to the TOF selection criteria and background subtraction 
scheme. The efficiency is close to 100% for the well-separated TOF peaks, but is only 
42% for the worst cases where the TOF peaks overlap. The contributions from inelastic 
scattering leading to the low excited levels for ^^^Cs and ^^''I are expected to be less than 
0.1% relative to the elastic processes. 

After all the correction factors are taken into account, the measured angular distri- 
bution for neutron elastic scatterings off Cs and I is displayed in Figure 0. The shape 
of the solid curve is the evaluated neutron elastic scattering cross sections on ^^^Cs and 
^^^I from ENDF/B-VI library |jT^ which comes from the Optical Model calculation. The 
absolute normalization is from a best-fit value. There is excellent agreement between the 



relative recoil angular distribution from this measurement and the predictions from the 
Optical Model. In particular, this is a kinematical regime where the size of the nuclei 
is comparable to the wavelength of the incident neutrons, and the diffraction pattern is 
well- reproduced by the data. 

The measurements demonstrate the validity of the nuclear recoil event selection pro- 
cedures at all angles. In addition, the results in Figure ^ confirms the Optical Model 
predictions on neutron elastic scatterings from a direct measurement of nuclear recoils 
on heavy nuclei (instead of from the scattered final-state neutrons in transmission-type 
experiments) — that is, they show that neutrons do elastically scatter off from the nuclei 
but not via some anomalous processes. 



4 Pulse Shape Discrimination 

The nuclear recoil measurements provide in addition an excellent data set for the study of 
PSD at small photo-electron level near detection threshold. Depicted in Figure ^ are the 
measured nuclear recoils at 50° and 95°, as well as the 7-pulse due to the 22.1 keV line from 
^°^Cd, averaged over many events and normalized to the same light coUecition. It can be 
seen that heavily-ionizing events such as those due to a-particles and nuclear recoils have 



faster decays than those from 7's — opposite to the response in liquid scintillator ||T2 . 

The neutron induced nuclear recoil pulse shapes are similar between 50° and 95° 
scattering angles, corresponding to 4.9 keV and 13.4 keV electron-equivalence light yield. 



respectively. The shape is also similar to that of the a source data |T^ taken at high 



energy. At smaller scattering angles, the accidental low-energy 7-background becomes 
important such that the intrinsic nuclear recoil pulse shape cannot be derived. In the 
PSD studies, the recoil pulse shapes for these data set were taken to be the same as those 
measured at 50°. The similarity of the pulse shape between 50° and 50° justifies this 
assumption. 

Two independent approaches to obtain PSD were studied: (1) neural network tech- 
nique [0 and (2) the "classical" method based on mean-time plus partial charge. For 
the neural network approach, we used 4000 events each from the above two data sets as 
training samples and additional 4000 events each as testing samples. The trained neural 
net has 20 input nodes (that is, 5 FADC time bins combined to form one node in Fig- 
ure ^ in 5 //s window with 10 hidden nodes. The input node value should be greater 
than zero after pedestal correction. A negative value is reset to zero and the sum of the 
20 input values are normalized to one. The algorithm used is back propagation. The 
excellent separation power from neural net study is depicted in Figure ^ The trained 
neural net is apnlied to the complete 50° data set. Good separation is obtained between 



the 43 keV nuclear recoils and the low energy 7 pulses. A cut on neural output greater 
than 0.9 selects almost all the Csl nuclear recoil events. The number of events as well 
as the mean and RMS electron-equivalence energy are consistent with those determined 
from the statistical background subtraction method in Section The PSD separation 
with neural network method does not apply to the 20° data set. The neutron-beam en- 
vironment gave rise to large accidental background (mostly low energy 7-events) which 
produce much stronger bias for the low energy nuclear recoil signatures. This would not 
be a problem in a realistic Dark Matter experiment where the background rate would be 
extremely low. 

A large simulated sample was generated at different photon-statistics situation with 
these recoil/7 pulse shapes, and the results are displayed in Figure ||. It can be seen 
that a 95% 7 rejection with >90% efficiency for pulses with an average 20 photo-electrons 
can be achieved. It is projected that the performance can be further improved in a low- 
background quiet environment which allows a longer sampling window. These are better 
than the PSD capabilities for Nal(Tl) detectors at the same light yield 0. 

The alternative and complementary approach is one based on mean-time and partial 
charge. Two variables, the mean time (t) and the partial charge ratio R, are defined for 
each pulse, such that 

50 

E (Aiti) 

(t) = ^^^50 (2) 

E Ai 

i=l 

and 

30 

E Ai 

R = '-^ — , (3) 
E Ai 

i=l 

where Ai is the FADC amplitude at time bin ti. The scattered plots for (t) versus TOF 
as well as (t) versus R for the data set at 95° are shown in Figures pUjafcb, respectively. 
The electron-equivalence energy for Csl recoils is 13.4 keV at this scattering angle. There 
are clear separations between the Cs and I nuclear recoil events, represented by the box 
region, and the accidental 7-background for events triggered by the recoils of the wrapping 
materials (protons as well as other heavy nuclei). The cuts in (t) versus R space as defined 
in Figure |lO|b are subsequently applied to the other data set at lower recoil energies to 
select nuclear recoil events. 

The events due to the 22 keV 7-rays from a ^°^Cd source are displayed in Figure p!0|c. 
The good overlap between the box and the events justifies that the cuts are valid to 
select "in-time" events. However, the separation between nuclear recoils from the in- 
time 7-background is not satisfactory. While the classical method is sufficient for event 



identification in this measurement where the main background are the accidental low- 
energy 7's, it would be inadequate for a Dark Matter experiment where the 7-background 
also provide the trigger and are therefore in-time. In comparison with Figure the neural 
network method gives 95% rejection of the in-time 7's at a signal efficiency of better than 
90% at >20 photo-electrons (that is, larger than 50° in this data set). The neural network 
method, therefore, will give superior performance in a Dark Matter experiment compared 
to the classical method. 

The electron-equivalence light yield at different scattering angles for the nuclear recoil 
events selected by either of the two PSD algorithms are measured. The comparison 
between the two PSD algorithms with the subtraction method discussed in Section |^ is 
depicted in Figure |Tl|, indicating that consistent results are obtained. The deviations 
of the mean values among these independent measurements are also consistent with the 
estimated systematic uncertainties. 

Therefore, the robustness of the quenching factor results in Figures ||a&b measured 
by the statistical subtraction method discussed in Section ^ is further established. This 
method involves a straight-forward normalization and subtraction scheme, and is therefore 
expected to give more unbiased results compared to the two PSD algorithms. 



5 Discussion and Summary 

In this article, we presented a new measurement on the quenching factors of Cs and I in a 
CsI(Tl) crystal scintillator based on statistical background subtraction. Lower threshold 
and improved accuracies are achieved compared to previous measurements. The measured 
differential cross-section of neutron scattering from I and Cs is in excellent agreement with 
Optical Model derivations, and represents the first confirmation of the Optical Model on 
neutron elastic scattering cross-sections with a direct measurement of nuclear recoils from 
heavy nuclei. The recoil differential cross sections and the quenching factors are relevant 
to the studies of radiation damage in materials. 

Two complementary pulse shaping discrimination techniques at the low photo-statistics 
regime were studied based on the nuclear recoil data. It is shown that the neural net- 
work approach gives superior rejection power to in-time 7-background compared to the 
classical method based on mean-time and partial charge. Both procedures are adequate 
for differentiating accidental 7's. The electron-equivalence light yield derived from both 
PSD methods are consistent with those obtained by statistical subtraction. The internal 
consistencies of the light yield measurements among the three different algorithms and 
also with the Optical Model predictions of differential cross section enhance the reliability 
and robustness of the results. 



By optimizing the detector geometry and using green-extended photo-cathode, a de- 
tector with several kg modular mass and a light yield of a few photo-electron per keV 
can be realized The background from ambient radioactivity or intrinsic radio-purity 
is of course crucial to all low-background experiments. Their considerations for CsI(Tl) 
detectors are discussed in Refs. 0, |2^. Levels of better than the 10~^^ g/g level in 
concentration for the ^ss-y ^j^^^ '^^'^Th series have been demonstrated, assuming secular 
equilibrium. The potential problem of the internal ^^''Cs contaminations can be over- 
come via selection of clear ore materials and careful chemical processing and purification 
treatment [|To|l. 



Displayed in Figure |T2| is the sensitivity plot for spin-independent WIMP interaction 
cross section per nucleon with a target mass of 1 ton of CsI(Tl). The sensitivities are 
based on the assumptions that the detection threshold is at 2 keV (or about 16 keV 
recoil energy) and still allows enough photo-electron statistics for background rejection 
with PSD, and that the observed background rate after applying the PSD cuts is at the 
same as or ten times better level than the Nal(Tl) experiments (1 and 0.1 per day per 
kg per keV, respectively). It can be seen that the DAMA allowed region, represented by 
the shaded region in Figure 0, can be thoroughly probed. Given that there is matured 
experience of scaling up the CsI(Tl) detector to multi-ton systems, there are potentials of 
further improvements in the sensitivities, and in the case of positive results, performing 
an accurate measurement of the annual modulation. 
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Table 1: Measured results of the neutron elastic scatterings from CsI(Tl) detector. Errors 
shown are combined systematic and statistical uncertainties. 




Figure 1: Schematic diagram of the experimental set-up. 




Figure 2: The TOF spectra at 50° scattering angle for (a) with wrapping material only, 
(b) with CsI(Tl) target. 




Figure 3: The energy spectra at 50° scattering angle for (a)events with high TOF values, 
(b) events with low TOF values. 




Figure 4: The observed energy spectra with background subtracted at different scattering 
angles. 
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Figure 5: The normalized pulse shape in the first 100 FADC samplings at 50 ns per bin. 
Solid historgram is due to 22.1 keV X-rays, while the two dashed histograms are from 
recoils at 50° and 95°. 




Figure 6: (a) The quenching factors and (b) electron-equivalence light yield versus recoil 
energy measured in this work shown as black circles. The solid line is a linear fit to the 
data. Open triangles and open squares are data from Ref. M and Ref. fl, respectively. 
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Figure 7: The differential cross sections determined from nuclear recoils in the target. 
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Figure 8: The neural network output of 4000 nuclear recoil events (solid histogram) and 
4000 7 events (dashed histogram). The nuclear recoil energy is 43 keV and the 7 energy 
is 22 keV. 




Figure 9: The PSD performance based on neural net studies. The 7 rejection is fixed 
at 95% and the result of signal efficiency is plotted against available number of photo- 
electrons. 
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Figure 10: Scattered plots for the (t) parameter versus the (a) TOF values and (b) R 
parameters for the nuclear recoil events at 95°, and (c) for events due to the 22 keV 7-rays 
from a ^°^Cd source. The box region represents the Csl nuclear recoil events, indicating 
clear seperations from the accidental 7-background. 
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Figure 11: Ratio of electron-equivalence light yield derived from the two PSD algorithms, 
labelled "Mean Time" and "Neural Net", with the measurements of the statistical sub- 
traction method. Consistent results are obtained among them. 
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Figure 12: The sensitivity plot of cross-section versus WIMP mass on a 1 ton CsI(Tl) 
detector, with the assumptions that the detector threshold is set at 2 keV with a quenching 
factor of 0.13. The two cases of background counting rate at 1 and 0.1 counts per day 
(cpd) per kg per keV are shown. The DAMA allowed region is presented by the shaded 
area. 



